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ABSTRACT 

We have used a time series of aerosol surface based on 
the measurements of Hofmann (1990) to investigate the 
modulation of total column ozone caused by the perturba- 
tion to gets phase chemistry by the reaction A^ 20 s(gas) -f 
H 2 0(euevo) — ► 2H N Oa(gas) on the surface of stratospheric 
aerosols. We have tested a range of values for its reaction 
probability, 7 — 0.02, 0.13 and 0.26 which we compared to 
unperturbed homogeneous chemistry. Our analysis spans 
a period from Jan. 1974 to Oct. 1994. The results sug- 
gest that if lower values of 7 are the norm then we would 
expect larger ozone losses for highly enhanced aerosol con- 
tent than for larger values of 7 . The ozone layer is more 
sensitive to the magnitude of the reaction probability un- 
der background conditions than during volcanically active 
periods. For most conditions, the conversion of NO 2 to 
HNO3 is saturated for reaction probability in the range of 
laboratory measurements, but is only absolutely saturated 
following major volcanic eruptions when the heterogeneous 
loss dominates the losses of N 2 Os> The ozone loss due to 
this heterogeneous reaction increases with the increasing 
chlorine load. Total ozone losses calculated are compara- 
ble to ozone losses reported from TOMS and Dobson data. 

The appearance of massive ozone depletion over 
Antarctica has prompted a great interest in heterogeneous 
chemistry on polar stratospheric clouds (PSCs). More 
recently, it has become apparent that heterogeneous re- 
actions may also be playing an important role in ozone 
destruction outside of the polar regions (Hofmann and 
Solomon, 1989; Brasseur et ah, 1990; Pitari et al., 1991; 
Rodriguez et ah, 1991). The ozone destruction from het- 
erogeneous reactions is related to the availability of aerosol 
surface which is modulated by potent volcanic eruptions 
which penetrate the stratosphere. 

In this study, we evaluate this modulation of ozone by 
volcanic eruptions using a one-dimensional photochemical 
model and input the aerosol surface as a time series based 
on the measurements presented by Hofmann (1990). We 
have assumed that the aerosols are spherical and can be 


represented by a log-normal size distribution. The mea- 
surements spanned a period from 1974 to April, 1990. 
Prior to 1974 and from the end of the data set to June 
15, 1991, the aerosol surface was kept constant to represent 
aerosol background conditions. To represent the impact of 
Pinatubo for the period from June 15, 1991 to October 13, 
1994, the aerosol parameters from the El Chichon pulse, 
taken to span the period from April 4, 1982 to April 1 st, 
1990, were used and the amount was scaled by 1.5 (Lar- 
son et al., 1982). Consequently, we have extended the time 
series from 1974 to late 1994 (Figure 1 ). The profiles of 
aerosol surface adopted were based on measured profiles 
of aerosol mass mixing ratio taken in mid- 1989 (Hofmann, 
1990). In our model, the aerosol layer extends from 7 to 
27 km with peak mixing ratio at 22 km and the profile is 
assumed to remain unchanged with time. 

They are a number of heterogeneous reactions of pos- 
sible importance. In this note, we concentrate on one het- 
erogeneous reaction: 

N 2 0$ (gas) + H 2 0( aero) - 2 HNO z (gas) ( 1 ) 

(Evans et al., 1985; Hofmann and Solomon, 1989). This 
reaction affects the partitioning of nitrogen species. It 
channels odd nitrogen from a short-lived reservoir species, 
N 2 Os, to a longer-lived one, HNO 3 . Although the de- 
struction of ozone by N0 2 is expected to decrease, it is 
expected to be ofFset by an increase in active Cl x (eg. Ro- 
driguez et al., 1991). 

We have used a simple parametrization for the het- 
erogeneous chemistry based on the aerosol surface area 
and measured reaction probability. Laboratory measure- 
ments of reaction probability ( 7 ) of A r 2 0 5 + fl 2 0 on 

aerosol surface range from 0.04 at 283° K (Van Doren et al., 
1990) to 0.09 and 0.139 at 293°K and 274°K respectively 
(Mozurkiewich and Calvert, 1988) and to 0.14 at 215°K 
(Hanson and Ravifehankara, 1991). These rates were found 
to be independent of the relative humidity (Mozurkiewich 
and Calvert, 1988). To account for possible temperature 
dependencies in the reaction probability and for uncertain- 
ties in total aerosol surface (emerging from assumptions 
made in its determination from published data), we have 
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Figure 1: Aerosol surface time series at aerosol peak mix- 
ing ratio. See text for details. The peaks represent vol- 
canic eruptions: 1. Fuego, 2. Mount St. Helens, 3. Alaid, 
4. Nyamuragira, 5. El Chichon, 6. Nevado del Ruiz, and 
7. Pinatubo. The Pinatubo pulse is assumed to be 1.5 x 
El Chichon’s (Larson et ah, 1992). 

chosen the reaction probability to cover a range of values: 
0.02, 0.13 and 0.26. We have compared the results of these 
runs to a homogeneous chemistry case. 

Our analysis spanned a period from January 1974 to 
October 1994. Scenarios were run with the chlorofluoro- 
carbons and halons fluxes increasing according to OECD 
estimates to 1980 (OECD, 1981), according to global pro- 
duction rates from 1980 to 1985 and follow both the 
Montreal Protocol and the London amendments thereafter 
(eg. McConnell et al., 1990). 

Figure 2 shows time series of total column ozone for 
various values of 7. For the period from early 1974 to 
late 1994, the loss of ozone due to homogeneous chemistry 
increases monotonically in response to increasing chlorine 
loading. Calculations for heterogeneous cases clearly show 
signatures of major volcanic eruptions especially those of 
Fuego (1974), El Chichon (1982) and Pinatubo (1991). 
Ozone loss maximizes for El Chichon ( — 16 to —19 DU) 
and Pinatubo ( — 14 to —19 DU) as compared to the total 
column ozone in January 1979 and January 1990, respec- 
tively. Maximum loss of ozone lags the maximum aerosol 
surface peak by 2 to 4 months. 

The effect of the heterogeneous reaction of N 2 0$ is 
determined by how close the chemistry is from satura- 
tion . By saturation it is meant that the rate of reac- 
tion N 2 Os (gas) + H 2 0 (aero) — ► 2 HNO 3 (gats) 

is limited by the production rate of A 2 0 s which is, in 
turn, determined by the nightime formation of A r 0 3 via 
NO 2 + O3 —► A r 0 3 + O 2 (with rate k 2 ). It is expected that 
most of the A 2 05 formed will be converted to //A r 0 3 , gen- 
erating a significant reduction in A r 0 2 densities and an 


associated increase in CIO concentration. This process 
increases the CIO concentration above that calculated for 
unperturbed homogeneous chemistry. The CIO amplifica- 
tion is limited by the amounts of available NO x which is 
determined by an equilibrium between the loss of NO x to 
N 2 0& formation and its production from HNO 3 photoly- 
sis. 

Figure 3 shows time series for the loss time constants 
for N 2 Os to form HNO 3 (ri) (with 7 = 0.13), to photolyse 
(r 2 = l/J/v 2 o 6 ), and for N0 2 loss to 0 3 (t 3 = l/fc 2 (0 3 )). 
For 7 = 0.13, the chemistry is saturated with N0 2 *f 
03 — ► NO 3 4- 0 2 , the slowest reaction, being rate limit- 
ing. It can be seen that, during periods of volcanic qui- 
escence (mid-1976 to 1980; 1988-91), photolysis of jV 2 0 5 
is more rapid than the heterogeneous loss of A r 2 0 6 , even 
for 7 = 0.13. During volcanically active periods (1980-86; 
1992-thereafter), the 7V 2 0 5 conversion to HNO 3 is faster 
than its photolysis. The size of 77 versus the length of the 
night ( rjv) is indicative of the strength of the NO 2 con- 
version to HNO 3 . If T\ < Tjv, most of the N 2 0$ formed 
during the night will be converted to HNO 3 and the nigh- 
time N0 2 loss will be maximized. On the other hand, if 
T\ > T/v, a fraction of the N 2 Os formed during the night 
will be present at sunrise. The heterogeneous reaction will 
then be competing with N 2 Os photolysis. For background 
aerosol conditions, photolysis dominates the N 2 Os losses 
and a large fraction of the N 2 0$ will be converted back 
to N0 2 . This limits the increase in 0/0 concentration 
and reduces the ozone loss. The magnitude of the sunrise 
concentration of ZV 2 0 5 is inversely proportional to the re- 
action probability making the chemistry more sensitive to 
7 during periods of volcanic quiescence (Figure 2). 

During volcanically active periods, 7 < 7 . We call 
this situation absolute saturation because by sunrise, 
JV 2 0 5 will have been mostly converted to HNO 3 with only 
a small fraction left to regenerate the N0 2 through the rel- 
atively slower A r 2 0g photolysis. Under these conditions, 
lower values of ozone column are found. The balance be- 
tween N 2 0$ H" /f 2 0( aeroso I) and N 2 0$ photolysis is deter- 
mined by the aerosol loading. Periods of enhanced aerosol 
content tend to favour the heterogeneous reaction leading 
to a less efficient, but still active, photolysis. For instance, 
in 1979, our calculations show that 64 % of the daily aver- 
aged N 2 Os was converted back to N0 2 by photolysis and 
only 34 % was channelled to HNO 3 by the heterogeneous 
reaction. At the peak of the El Chichon aerosol pertur- 
bation (1983), only 4 % of the N 2 0 5 was converted back 
to N0 2 while 96 % was transformed into HNO 3 . What 
emerges from this work is that even for saturated chem- 
istry the aerosol loading will modulate the N0 2 and CIO 
concentrations by controlling the relative balance between 
heterogeneous conversion of N 2 0*> and its photolysis. This 
can be seen in the time series of total column of ozone (Fig- 
ure 2) as the well defined troughs matching the aerosol 
peaks of Figure 1. It can further be added that for en- 
hanced aerosol levels, the chemistry becomes less sensitive 
to changes in aerosol loading. 
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Figure 2: Total column ozone time series from 1974 to 1995. Cases include homogeneous 
chemistry and the heterogeneous reaction N 2 0$ (gas) + H 2 0 (aero) — + 2HN0$ (gas) 
with 7 = 0.02, 0.13, 0.26. 


The magnitude of the amplification of CIO densities 
above the background levels determines how much ozone 
will be lost with the inclusion of reaction (1): the larger 
the reaction probability, the lower the possible CIO am- 
plification following volcanic eruptions. For instance, in 
the lower range of values, for 7 = 0.02, the ozone is more 
sensitive to aerosol loading since there are no, or at least re- 
duced, saturation effects. For this case, total column ozone 
changes were calculated to be as large as — 5.8% between 
1979 (aerosol background conditions) and 1983 (El Chi- 
chon’s maximum ozone loss). For the same period changes 
of —5.3% and —5.0% were calculated for 7 = 0.13 and 
7 = 0.26 respectively. For the Pinatubo eruption, total 
ozone changes between June 1990 and March 1992 were 
calculated to be: —6.0%, —5.2% and —4.6% for 7 = 0.02, 
0.13, 0.26, respectively. For the larger values of 7, abso- 
lute saturation will set in sooner, the time depending on 
the aerosol surface area, and the apparent ozone loss will 
not be as great as shown by our calculations. 

Following volcanic eruptions, the temporal behaviour 
of the recovery curve is slower for the higher 7 values. For 
the low 7 values, the recovery is more rapid and the ef- 
fects of the diminished levels of aerosols are seen in ozone 
recovery with a short time lag. For the large values, the re- 
covery is much slower until the aerosol levels become small 
enough to lower the efficiency of the heterogeneous reac- 
tion beyond the level of absolute saturation. For example, 
the curves with 7 = 0.02 have a 50 % recovery by 1.5 
years while it takes about 3 years for the case of 7 = 0.26. 
Our calculations show that for reaction probability 


within the range of measured values, the ozone loss be- 
tween January 1979 and May 1990 is < 3.5% (Table 1). 
These values are consistent the 4 % average loss observed 
by TOMS in the 20-60° N area (Stolarski et ah, 1992; 
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Figure 3: Time series of loss time constants at 20 km 
from 1974 to Oct. 1994 for NO 2 lost through the reac- 
tion N0 2 + 0 3 — + NO 3 4- 0 2 ; N 2 O b lost by N 2 0 5 (gas) + 
H 2 0(aero) — ► 2//A r 03(gas) with 7 = 0.13; and N 2 O s lost 
by photolysis. 
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7 

0.00 

0.02 

0.13 

0.26 

O 3 loss (%) 

-2.1 

-2.4 

-3.1 

-3.5 


Table 1: Total column ozone loss from Jan. 1979 to Jan. 
1989. 


7 

0.02 

0.13 

0.26 

CIO amplification 

1.3 

2.1 

2.5 


Table 2: CIO amplification factor vs 7 The amplification 
factors were calculated at 20 km as the ratio of the CIO 
concentrations for the heterogeneous to the unperturbed 
homogeneous cases. For given 7, the CIO amplification 
factors calculated in Jan. 1979 and Jan. 1990 are identical. 

Stolarski et al., 1991). It is also within the limits of the 4 
% loss from 1980 to 1990 reported for the North- American 
stations of the Dobson network (Stolarski et al., 1992). 

Anthropogenic emissions of CFC and halon compounds 
enhance the ozone loss due to reaction (1). For instance, 
if 1979 and 1990 are assumed to be background years for 
aerosol loadings and if we use January 1948 as a reference 
year (334 DU), our calculations show that the ozone loss 
due to reaction (1) with a 7 = 0.13 was 2.4 % in 1979 
and 3.4 % in 1990. The enhanced ozone loss is due mainly 
to an increase in atmospheric chlorine load. It was found 
that for lower aerosol levels such as is observed during 
aerosol background years, the CIO concentration is en- 
hanced above the unperturbed homogeneous chemistry by 
a fixed factor which depends on the reaction probability of 
reaction (1) (see Table 2). For instance, for 7 = 0.13, the 
CIO concentrations in 1979 and in 1990 were 2.1 times the 
CIO concentrations calculated for the homogeneous chem- 
istry case. Thus, the proportion of the ozone loss will 
increase as the atmospheric chlorine load increases. Con- 
sequently, under aerosol background conditions, a large 
reaction probability will generate a larger departure from 
the unperturbed homogeneous chemistry case. 

Compared to unperturbed homogeneous chemistry, the 
heterogeneous reaction A^C^T^C^aero) — ► 2//AfOa(gas) 
lowers the total column of ozone. Following large volcanic 
eruptions (eg El Chichon and Pinatubo), maximum losses 
of 19 DU were calculated with maximum relative losses 
close to 6 % for the lower reaction probability. Recov- 
ery from these eruptions is longer for the larger reaction 
probability. During background aerosol years, the ozone 
layer appears more sensitive to the magnitude of the re- 
action probability with calculated changes in total column 
of ozone proportional to the reaction probability. These 
background years correspond to a chemistry which is not 
absolutely saturated, leaving room for extra NO x conver- 
sion to HNO3. The 1979-1990 trends in total column of 
ozone show a small departure from the unperturbed chem- 
istry and remain within the range of observed TOMS 


and Dobson data. Finally, the impact on the oz<5ne layer of 
the reaction studied amplifies with the increasing chlorine 
loadings. 
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